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We study the effect of rotation on global properties of neutron star with a hyperon core in 
an effective chiral model with varying nucleon effective mass within a mean field approach. The 
resulting gross properties of the rotating compact star sequences are then compared and analyzed 
with other theoretical predictions and observations from neutron stars. The maximum mass of the 
compact star predicted by the model lies in the range (1.4 — 2.4) Mq at Kepler frequency ^k, 
which is consistent with recent observation of high mass stars thereby reflecting the sensitivity of 
the underlying nucleon effective mass in the dense matter EoS. We also discuss the implications of 
the experimental constraints from the flow data from heavy-ion collisions on the global properties 
of the rotating neutron stars. 

PACS numbers: 21.65.+f, 13.75.Cs, 97.60.Jd, 21.30.Fe, 25.75.-q, 26.60.+C 



I. INTRODUCTION 

Recent advances in the study of cold dense matter has 
received new impetus both from observational astronomy 
and the recent heavy-ion collision experiment as well as 
the upcoming compressed baryonic matter (CBM) exper- 
iments, that can put constraints on the equation of state 
(EoS) of dense matter Arguably the uncertainties in 
the dense matter EoS can also be reduced appreciably 
through theoretical modeling of neutron stars and by an- 
alyzing their global properties with an eye on the obser- 
vational data from these compact objects. Consequently, 
building an EoS in accordance with the heavy-ion col- 
lision data and its subsequent application to the dense 
matter problem seems to be prerequisite. The study can 
throw significant light on potential problems such as the 
isospin dependence of nuclear forces, the presence of sta- 
ble hyperon in neutron star core and its impact on cold 
dense matter 0, 0] ■ 

Lately, there has been some striking signatures or ob- 
servations of both massive and the low mass stars, which 
may have interesting implications on the structure and 
composition of neutron stars [^]. For example, the ob- 
servation of massive compact stars (M « 2M0) 0, Q 
from observation of QPO's (Quasi Periodic Oscillations) 
from X-ray emmisions on one hand and also the least 
massive ones 1.18 ± O.O2M0 0] from the binary pulsar 
J1756 — 2251 emphasize the need to address key issues 
in limiting the observed global properties of the stars in 
t\i M — R plane. On the contrary, the canonical value 
of M = I.44M0, the largest precisely measured mass of 
the Neutron star PSRlQli -I- 16 Q effectively doesn't 
seem to provide any stringent condition and most of the 
existing realistic equation of state satisfies the criteria. 
These aforesaid facts signifies the need to constraint the 
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nuclear EoS of dense matter which are extrapolations of 
the properties defined in the vicinity of nuclear satura- 
tion density and therefore it is important to test these 
extrapolations with regard to the available experimental 
flow data from heavy-ion collisions. 

Equivalently considerable emphasis is being laid to 
understand nuclear matter EoS in the density range of 
(1 — 3)poj the region from where the entire possibilities of 
exotic forms of matter (hyperons, condensates of bosons 
or quarks) Il3l | starts appearing. It is known 

that the neutron star radius is sensitive to the properties 
of nuclear force in the immediate vicinity of nuclear sat- 
uration density (po)i whereas the maximum mass of the 
star is regulated by the dense matter EoS well beyond 
2po Among the various scenarios related to the star 
structure, the presence of stable hyperons is an exciting 
possibility and is known to have significant impact on 
the dense matter equation of state and in turn gets man- 
ifested in the star properties too [13, [3j The hy- 
perons such as A0(1116), i;-^0'+(1193) and S-'0(1318) 
appear at the expense of nucleons and forms a sizeable 
population in the density region of current interest. Al- 
though hypernuclei experiments [l6| supports the pres- 
ence of hyperons in dense matter, however this domain 
is also marred by uncertainties revolving around the na- 
ture of hyperon-hyperon interactions and their coupling 
strengths. Since they are the dominant species in the 
dense matter population, there is even greater need to 
understand the nature and behavior of these interactions. 

Apart from that, it is also known that the addition of 
quartic terms softens the EoS considerably in addition 
to that caused by the hyperons, which thereby makes 
it difficult to obtain a neutron star mass larger than 
2M0 (l7| . Here in the present investigation, we employ 
a model which has chiral symmetry and also embodies 
a dynamical generation of the vector meson mass along 
with the non-linear term in the scalar field interaction. 
The model has been used extensively to dense matter 
studies [isl . [l9l . [20| . which was then modified with the 
addition of non-linear terms in the scalar field interac- 
tion to reproduce nuclear saturation properties at rea- 
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sonable incompressibility (21I . [23 | . Finally we generalized 
the model to include the lowest lying octet of baryons in 
our previous work psj and applied it to study neutron 
star matter with varying incompressibility and var yin g 
nucleon effective mass. From our previous analysis [23|, 
it was apparent that the difference in the nuclear incom- 
pressibility was neither prominent in the EoS, nor in the 
static star properties, whereas the EoS was extremely 
sensitive to the underlying nucleon effective masses and 
hence was also visible in the global properties of the star. 
In continuation of our earlier work, we now extend our 
analysis to study the rotational attributes of neutron star 
with varying nucleon effective mass in the mean-field ap- 
proach within the same hadronic model. In short, our 
motivation here has been to investigate how the rota- 
tional attributes of the neutron star restrict the parame- 
ters of the present model. 

The outline of the paper is as follows: In the next sec- 
tion, we give a brief description of the ingredients of the 
hadronic model that we implement in our present work. 
After describing the numerical scheme for calculating the 
rotating neutron star sequences, the general features of 



the equation of state of hyperon rich matter is studied 
and the global properties of the rotating neutron star is 
presented. We shall analyze these results with respect 
to varying nucleon effective mass and also talk of a few 
constraints on the neutron star mass and radius imposed 
by recent observations. Finally we conclude with outlook 
on the possible extensions of the current approach. 



II. THE EQUATION OF STATE 

In this context, we tried to explore the consequences 
of an EoS based on chiral a model with nonlinear in- 
teraction terms in the scalar field. To have a realistic 
description of the dense neutron star matter, we con- 
sider the effective Lagrangian of the chiral model gen- 
eralized to include the lowest lying octet of baryons 
(ri,p, A°, 1]^'°'+, interacting through the exchange 

of the pseudo-scalar meson tt, the scalar meson a, the vec- 
tor meson uj and the iso- vector p— meson, and is given by 
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The first line of the above Lagrangian represents the 
interaction of baryons ^'s with the aforesaid mesons. In 
the second line we have the kinetic and the non-linear 
terms in the pseudoscalar-isovector pion field 'tt', the 
scalar field 'cr', and with x'^ = w'^ + a^. Finally in 
the last line, we have the field strength and the mass 
term for the vector field 'w' and the iso-vector field 'fP 
meson. The terms in eqn. (1) with the subscript 'B' 
should be interpreted as sum over the states of all bary- 
onic octets. In this paper we shall be concerned only 
with the normal non-pion condensed state of matter, so 
we take < if >— 0. 

The interaction of the scalar and the pseudoscalar 
mesons with the vector boson generates a dynamical mass 
for the vector bosons through spontaneous breaking of 
the chiral symmetry with scalar field getting the vacuum 
expectation value xq. Then the masses of the baryons, 
the scalar and the vector mesons, are respectively given 
by 



rriB = gaBXo, m„ = V2Xxo, = gu^BXo . (2) 



In the above, xq is the vacuum expectation value of 
the (J field. We could have taken an interaction of the 
p— meson with the scalar and the pseudoscalar mesons 
similar to the omega meson. However, a dynamical 
mass generation mechanism of the p— meson in a simi- 
lar manner will not generate the correct symmetry en- 
ergy. Therefore, we have taken an explicit mass term for 
the isovector p— meson similar to what was considered in 
earlier works [13, HU, HI, [H . 

We employ the mean-field procedure to evaluate the 
meson fields in our present calculations. In the mean- 
field treatment, one assumes the mesonic fields to be 
uniform i.e., without any quantum fluctuations. We re- 
call here that this approach has been extensively used 
to obtain field-theoretical EoS for high density matter 
p^ . and gets increasingly valid when the source terms 
are large [2J]. The details of the model that we use in 
our present investigation and its attributes such as the 
derivation of the equation of motion of the meson fields 
and its equation of state (e & P) of the many baryonic 
system, can be found in our preceding work [23| . For the 
sake of completeness however we write down the meson 
field equations. The vector and the iso-vector fields are 
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respectively given by 
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PB 



(3) 



In the above equations the quantity ps is the Baryon 
density and I^b is the 3rd-coniponent of the isospin of 
each baryon species. 
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The scalar field equation can be written in terms of the 



(4) variable Y = x/xq with x — {< a + n >) 
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where the effective mass of the baryonic species is = 
YruB and CaB = qIb/'^I are the c^jb = 5wb/"^w are the 
usual scalar and vector coupling constants respectively. 
It may be noted that the parameter 'A' in the Lagrangian 
does not appear explicitly in eqn. (5), but enters implic- 
itly through the mass term for the scalar meson, following 
equation (2). Similarly, in the present model describing 
dense matter, the meson mass is generated dynami- 
cally. This vector meson mass enters in Eq. ^ through 
the ratio c^j = [gui/mSY = l/a^o^ which has to be fixed 
alongwith Cc so as to satisfy the nuclear matter satura- 
tion properties, similar to what has been done in earlier 
works 0, [m, [23, ill. Further in Eq.®, the quantities 
Pb and psB are the baryon density and the scalar density 
for a given baryon species given respectively as. 
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where fc^ is the fermi momentum of the baryon and 7 = 2 
is the spin degeneracy factor. 

We now go directly to the total energy density 'e' and 
pressure 'P' for a given baryon density in terms of the 
dimensionless variable Y — x/x^ which is given as: 
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The terms in eqns. (3) and (4) with the subscript '_B' of total baryon number and charge neutrality condition 
should be interpreted as sum over all the states of the given by [23| 
baryonic octets. The meson field equations for the a, uj 
and p— mesons are then solved self-consistently at a fixed 
baryon density to obtain the respective field strengths. 
The EoS for the /3— equilibrated for the hyperon rich 

matter is obtained with the requirements of conservation /-^ QbPb + QiPi — 0, (10) 

B I 
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where pB and pi are the baryon and the lepton (e,/i) num- 
ber densities with Qg and Qi as their respective electric 
charges. 

Using the computed EoS for the neutron star se- 
quences, we calculate the structural properties of neutron 
stars with a hyperons core. 



III. STELLAR EQUATIONS 

The equations for the structure of a relativistic spher- 
ical and static star composed of a perfect fluid were de- 
rived from Einstein's equations by Tolman, Oppenheimer 
and Volkoff [251, which are 



dP 
dr 



G[e + P] [M + 47rr3p] 



dM 
dr 



{r - 2GM) 



(11) 



(12) 



with G as the gravitational constant and M(r) as the 
enclosed gravitational mass. We have used c = 1. Given 
an EOS, these equations can be integrated from the origin 
as an initial value problem for a given choice of central 
energy density, [sc)- The value of r (= R), where the 
pressure vanishes defines the surface of the star. We solve 
the above equations to study the structural properties 
of a static neutron star using the EoS derived for the 
electrically charge neutral hyperonic dense matter [26l . 

US- 

Although a relativistic compact star has much com- 
plicated internal structure but its properties can be rea- 
sonably approximated by some simplifying assumptions. 
The matter inside is assumed to be a perfect fiuid on 
the basis of the observation from pulsar glitches, which 
shows that the departure from perfect fluid equilibrium 
due to the solid crust is quite negligible (~ 10~^) [2^ . 
At birth, a neutron star is differentially rotating, but be- 
cause of several factors such as the cooling phenomenon, 
shear viscosity, neutrino diffusion etc., the star assumes 
a uniform rotation. So the approximation of a zero tem- 
perature perfect fluid neutron star matter is a good one. 

In order to calculate the models of rotating star basi- 
cally two approaches are being employed [29|, namely 
the Hartle approach (slow rotation) [3011 and the Ko- 
matsu, Eriguchi, and Hachisu (KEH) method (fast rota- 
tion) [3l|, [34I approach. Although an improved version of 
the former approach was em ploy ed to calculate the prop- 
erties of rotating stars [ssl |34 [ssj . large descripancies 
has been noted compared to models without the assump- 
tion of slow rotation, particularly near the mass-shedding 
limit 36]. In our present investigation, we employ the 
KEH approach to calculate the model for rapidly rotat- 
ing stars near the mass-shedding limit or conversely upto 
the Kepler frequency, which we briefly describe now. 

Let us consider equilibrium stars in uniform rota- 
tion with static, axial symmetric space-time. Now the 



time translational invariant and axial-rotational invari- 
ant metric in spherical polar coordinates {t,r, 9,(f>) can 
be written as 



ds^ = ~e-"'dt-'+e''"{dr^+r^d0^)+e'"'r^ sin-" e{d(^-Ludtf, 

(13) 

where the metric functions v,a,(3,oj depends only on r 
and 0. By treating stellar matter as perfect fluid, the 
energy momentum tensor can be written as 



with the four-velocity 



Vl — w' 



:(1, 0,0,17). 



Here 



(14) 



(15) 



(16) 



is the proper velocity relative to an observer with zero 
angular velocity and 17 is the angular velocity of the star 
measured from infinity. Now we can compute the Ein- 
stein field equations given by 



(17) 



(where R^y is the Ricci tensor and and R is the scalar 
curvature and with c = G = 1). This leads to the equa- 
tions for the metric functions as 



A- 



1 d 



A 

1 d 



r dr r"^^ dp 



pe2 
7e' 



A 



2_a 

r dr r^ 



d_ 

"dp 



where 



A 



dr^ 



2 d 



1 92 



r dr r^ 86^ 



de 



Spir,p), 


(18) 


S^{r,p), 
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(20) 
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7 — f3 + i'jP — V ~ and p = cos 0. We refer [3l| for 
the explicit expressions of the source terms Sp and Sui ■ 
These basic differential equations can be transformed into 
an integral representation to handle the boundary condi- 
tions in relatively simple manner with the help of multi 
dimensional Green's functions and solved by Komatsu- 
Eriguchi-Hachisu method[3ll. [3^. 

A maximum limit for the stable rotation of a star fix, 
is set by the onset of mass shedding from the equator 
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of the star. General relativistic expression for this Kep- 
lerian frequency '^Ik' can be obtained by using the ex- 
tremal principle to the circular orbit of a point particle 
rotating at the equator of the star [s?] . This leads to the 
expression for the Kepler frequency, which is given as, 



1 I' 



2i}' 



1/2 



, (21) 



where, ip — f3 + and ' denotes the differentiation 
with respect to the radial co-ordinate r. 

Next, the general relativistic expression for the mo- 
ment of inertia / of a rotating star is given by [ssj 



/ = ^ / drdOdcj^^-gTl, 



(22) 



where g is the determinant of the metric (|13p . Using 
eqn. (14-16), the moment of Inertia is then given as, 



I=?^JdrJ ddr^ sin' 9^-^^^ e'^'^+^\ (23) 

and hence the angular velocity 'J' of the star can be 
obtained from the relation 



J = m. 



(24) 



We use the code written by Stergioulas |39| based on 
the Komatsu-Eriguchi-Hachisu method to construct uni- 
formly rotating star models. 



IV. RESULTS AND DISCUSSIONS 

A. Choice of parameters and compatibility with 
the Heavy-Ion collision data 

From our earlier work [2^ it was conclusive that the 
effect of varying incompressibility on the underlying nu- 
clear equation of state and the star properties is mini- 
mal. There we found that EoS with incompressibilities 
K = 210, 300 & 380 MeV resulted in stars with nearly 
same maximum mass at similar central densities. Simi- 
lar results were obtained in other global properties of the 
compact star such as gravitational redshift, maximum 
baryonic mass and radius too. However the difference in 
nuclcon effective mass was found to be prominent both 
in the dense matter EoS and also on the star proper- 
ties. We now investigate the effect of varying nucleon 
effective mass on the rotational attributes of the evolved 
neutron star using the same model and same set of pa- 
rameters. It is to be noted that the parameter set is 
in accordance with recently obtained heavy-ion collision 



data [40| ■ Therefore in the present context it shall be in- 
teresting to see the implications of the constraints from 
heavy-ion collision data on neutron star structure. 

Here it is imperative to mention that we chose the pa- 
rameters for the present analysis keeping the incompress- 
ibility {K — SOOMeV) constant. The incompressibility is 
higher in comparison to other successful mean field mod- 
els like the NL3 {K = 21l.bMeV) [4l|. However, is well 
known that the nuclear symmetry energy greatly influ- 
ences the EoS at high densities. The value for symme- 
try energy J — 37AMeV in case of NL3 is unphysically 
large in comparison to that we have in the present model 
(J = 32MeV). On the other hand the DBHF, which is 
considered to be more realistic EoS in the non-relativistic 
domain has a symmetry energy of J = 31.5MeV and is 
comparable to our prescription and is in good agreement 
with the empirical models and also with the variational 
approach f43|. In addition to these aforesaid differences 
in the nuclear saturation properties, there are variations 
in the nucleon effective mass in various models. We take 
TO*/m„ = (0.80 — 0.90) which is consistent with the 
values obtained from the analysis of neutron scattering 
off lead nuclei IIH. 
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FIG. 1: Pressure {MeV fm,-^) as a function of normalized 
baryon density (ps/po) of symmetric nuclear matter. The 
shaded region is consistent with the experimental flow data 

m- 

The parameter set for the present work is listed 
in Table- 1, which vary in nucleon effective masses 
{m*jq/mN = 0.80,0.85 & 0.90), but has the same incom- 
pressibility value dX K = ZQQMeV . The parameter sets 
satisfies the nuclear saturation properties, Eb-, energy per 
nucleon, —16.3 MeV at saturation density 0.153 fm~^, 
effective nucleon Landau mass (0.8 — 0.9) totv, incom- 
pressibility K = 300Mey, and asymmetry energy coef- 
ficient value (~ 32 MeV), so that our extrapolation to 
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TABLE I: Parameter sets for the model. 
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FIG. 2: Pressure {MeV fm''^) as a function of normalized 
baryon density (pn/po) of neutron matter. The shaded re- 
gion is consistent with the experimental How data [43|. The 
upper polygon represents the stiff EoS, and the lower one cor- 
responds to data for a soft EoS. 



higher density remains meaningful. The coupHng con- 
stant CpN is fixed to reproduce agym to the empirical 
value, 32 ± 6 MeV[4|. This gives CpN = 4.66 fm^ 



for asyin=32 MeV. The first three columns in table 1 
gives the coupling strength of the meson fields and the 
nucleon, namely the scalar field c^n, the vector field 
Ci^AT and the iso- vector coupling strength CpN, where 
Cijv = {9iN /iTT-i)'^ , for i = <7,uj,p. The next two 
columns are the higher order scalar field constants 'B' 
and 'C, required to reproduce the bulk nuclear matter 
saturation properties with reasonable incompressibility 
which is given in the next column, followed by the nu- 
cleon effective mass for the three sets. With these pa- 
rameterization, TTio- = nipf^ {cuiN /co-n), turns out to be 
in the range (509 - 627)MeV for the three parameter sets 
enlisted in Table I. Note that in these calculations, other 
meson masses or the respective couplings do not enter 
separately but through the ratio's Cpjv — (ffpAr/wp)^ and 

It is to be noted that the nucleon effective mass taken 
in the present work lies in the range m^/mjv = (0.80 — 
0.90). These values are considerable larger as compared 
to the Walecka model. A higher effective nucleon mass 
is known to generate problems in describing finite nuclei 
properties. However, these values are consistent with the 
results derived from non-relativistic analysis of scatter- 
ing of neutrons from lead nuclei [43] . We shall show that 
the variation in nucleon effective mass shall have consid- 
erable impact on the underlying EoS of the dense matter 
regarding its composition and hence on the gross stellar 
properties. As stated earlier that the EoS considered for 
the present study is in accordance with recently obtained 
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heavy- ion collision data. In Fig. I, we compare our EoS 
with the experimental flow data from recent heavy-ion 
collision for symmetric nuclear matter [i^ . Here we find 
that among the three parameter set set II and III are in 
very good agreement within the entire density domain 
of « (2 — 5)po- Similarly in Fig. 2, we compare the 
case for neutron matter with the experimental flow data. 
The two shaded region correspond to the density depen- 
dence of symmetry energy. The upper one corresponds to 
a strong density dependence (stiff EoS) while the lower 
one corresponds to a weak density dependence (Soft EoS) 
of the symmetry energy ^40, J^]. Here, all the three EoS 
seems to be in good agreement with the experimental 
data. However set I seems represents the stiffest EoS 
among the three and set III is the softest prescription. 
Therefore, in the present context, it can be seen that the 
flow data favors a soft EoS. 



B. Composition of Hyperon rich matter 

The composition of charge neutral dense matter is 
very sensitive to the hyperon-meson coupling parameters 
which however are very poorly known [ISl . l46j . There are 
different phenomenological prescriptions e.g., using the 
quark counting arguments based on SU(6) theory to fix 
the vector couplings, while, fixing the scalar couplings 
from the potential depths of hyperon in nuclear matter 
(Tol [ssl . [431 . From various analysis the choice of the ratio 
of hyperon to nucleon coupling for a meson Xa < 0.72 has 
been emphasized IEJI and also from studies based on hy- 
pernuclear levels [4^, the choice (xq- < 0.9) is bounded 
from above. Following this convention, in our present 
work we take Xa^gcH / gcN^O.l , Xu,=gujH / 
and Xui —Xp, to calculate the EOS for the neutron star 
matter. This leads to the binding of A° in nuclear mat- 
ter: {B/A)t^=Xu,gu,^^Q + m*^ - ruA -30 MeV. We fur- 



ther assume that couplings to the E and S resonances 
are equal to those of the A hyperon [1, Q . It is worth 
noticing that taking the p— meson coupling as Xi^ =Xp 
or Xo- =Xp does not alter the EoS substantially. Further 
that among all the three coupling schemes, this choice 
restricts the equation of state of neutron star matter fol- 
lowing the constraint of A° binding in nuclear matter, 
ably supported by the hypernuclei experiments [l^. We 
would like to mention here however that the values of 
these potential are somewhat different [4^ from the anal- 
ysis of S~ atomic data and the final state interaction of 



S hyperons in experiments E22A at KEK as well as 
i?885 at AGS [SlJ, which indicates a repulsive potential 
for S~ hyperon in dense matter. However, the exact val- 
ues of the relativistic potential remain inconclusive. We 
shall carry out our numerical analysis with the values of 
the potentials as prescribed above. 

As stated earlier, the equation of state of dense matter 
is very sensitive to the nucleon effective mass. In Fig. 3, 
we show the respective particle densities of a charge neu- 
tral hyperon rich matter for m* /m = 0.80 (set I) upto ten 
times nuclear matter density. Similarly Fig. 4 and Fig. 5 
displays the particle population for set II {m* /m = 0.85) 
and set III {m* /m — 0.90) respectively. From the plots, 
the sensitivity of the nucleon effective mass on the re- 
spective particle composition is clearly reflected. It can 
be seen that with increase in the nucleon effective mass 
the appearance of various hyperons is pushed further to 
higher densities. The appearance or concentration of the 
hyperon species at a particular baryon density is known 
to regulate the degree of softening effect on the equation 
of state, which is then manifested in the global proper- 
ties of compact stars. For example, parameter set I and 
set II accommodates all the octet of baryons upto ten 
times normal nuclear matter density although the hyper- 
ons starts appearing at higher densities in the later case. 
However in case of set HI, the nucleon chemical poten- 
tial is not high enough to have E^" and 'EP even upto 
highest of density region. From the systematics of the 
nucleon mass in hyperon rich matter, we recall that set 
II is endowed with repulsive force dominance unlike set 
HI, which is dominated by the scalar forces [1^ . 
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FIG. 3: Particle densities (/m ^) in neutron star matter for 
parameter set I. 

In general, for all the parameter sets, the hyperons 
starts appearing at « 2po and are the dominant species 
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in par with nucleons at ~ lOpo- It can also be seen that 
the negatively charged particle species are favored in the 
dense matter composition followed by the neutral ones 
and the positively charged counterparts. For example in 
all the three cases, one finds that appears before A" 
although the vacuum mass difference between the two is 
« llMeV and the former is massive than the later. This 
is due to the fact that as the density is increased (« 2po), 
the negatively charged S~ starts competing with the neg- 
atively charged leptons in maintaining charge neutrality 
of the dense matter. Due to the depletion of the lep- 
ton concentration the corresponding chemical potential 
/ie decreases, thereby lowering the threshold chemical po- 
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tential of S as compared to A" . Similar arguments exist 
for the other negatively charged baryons. 



C. Global properties of a rotating Neutron star 



1 

/ 

/ 


I 
II 


1 1 1 1 




m 




''///'' 

H / 
('// 






1 




Rotating (o) = Q.^ 

Static (CO = 0) 

1,1 



I \ \ \ \ I I I L_ 

0.4 0.8 1.2 1.6 

£^ (lo'^g/cm^) 



FIG. 6: Dependence of Gravitational mass of the star (in Mq) 
on the central density of the star (lO^^grcm"^) for the three 
cases in comparison with the static results (solid lines). 

The properties of the neutron star is unique to the 
EoS considered. Using these three EoS we now calculate 
some of the global properties of the neutron star, which 
is rotating at its Kepler velocity (0_r-) by employing the 
RNS code [s^, as described in section III. 

In Fig. 6, we compare the variation of gravitational 
mass of the star for the three cases as a function of the 
central density of the star for both the static and the 
rotating one. In later case, we considered rotation till 
the onset of mass shedding at the equator or equiva- 
lently upto its Kepler frequency. The effect of rotation 
on the stellar mass is clearly visible from the plot, where 
it is noticeable that a rotating star acquires more mass 
than its static counterpart. For all the three case at 
(w = fiif), we notice an increase of nearly 17-20 % in 
the star mass in comparison to the static ones. However, 
this increase in mass doesn't alter the central density 
of the star very much. The maximum mass obtained 
for set I, II and III, for the rotating case are 2.4M0, 
1.9M0 and \AMq respectively and the central densi- 
ties obtained for both set I and II are 7.0 x lO^^^cm^'^ 
and for set III, it is 8.0 x lO^^'^cm^'^. Thus we find 
that for the three cases the maximum mass is obtained 
at a central density of (2.5 — 3)po- In case of non- 
rotating star the central density obtained at maximum 
mass lies in the range (3 — 3.5)po- The corresponding 



maximum baryonic masses for set I, II and III for the 
rotating case are 2.56Mq, 2.IIM0 and 1.47M0 respec- 
tively. The difference between the gravitational mass 
and the baryonic mass gives us the gravitational bind- 
ing of the star. The fact that at any given density the 
baryonic mass exceeds the gravitational mass of the star 
is a typical of compact stars. Consequently, we find 
that the EoS with lower nucleon effective mass results 
in more bound compact systems. From the observa- 
tional point of view, we find that recent observations of 
neutron star masses like M707 51+1807 =2.1±0.2 MqI^, 

2.0±0.1 M0[5^, MyeiaX-l = 1.86±0.16 

and Mye/aJf-2 = 1.78±0.23 [H predicts 

massive stars. Our results are in very good agreement 
with these observations, except for set III which results 
in low mass star among the three prescription and has the 
highest nucleon effective mass value. It is worth to re- 
call that although there are large observational errors in 
the mass-radius determination, yet in case of VelaX — 1 
the lower mass limit « (1.6 — 1.7)M0 is atleast mildly 
constrained by geometry (56j . 



2.5 



2 



' 1 ' 1 ' 1 ' 


1 






PSRJ0751 + 1807| ^\ 

■0-. \ 


E 

•V 






\ \ 






t \ 

^ Typical NS Mass •. \ - 
\t"--~ ■•• \ 


SET I 

- ■ ■ ■ ■ SET 11 

SETHI 

1 . 


1 







I I \ I I I I I \ I I I I 

22 24 26 28 30 32 34 
Radius (km) 



FIG. 7: Gravitational mass of the star (Mq) rotating at Ke- 
pler frequency as a function of the circumferential radius of 
the star (in km). Also plotted is the typical neutron star mass 
range (1.0 — 1.5) Mq and the mass estimates obtained from 
PSRJ0751 + 1807 with error estimate of la (2.1 ± 0.2) Me 
and 2a {2.lto.t)MQ ^. 

Of the three EoS, in the present analysis, parameter 
set I is the stiffest prescription and set III is the softest 
one. This property is well documented in Fig. 7, where 
we plot the mass of the neutron star as a function of 
the radius and compare our results with the typical ob- 
served mass of neutron stars (1.0— 1.5)M0 and also with 
more massive ones such as the case of PSR J0751 -I- 1807 
within the error estimates of la and 2a. Here we find 
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FIG. 8: Moment of Inertia {W'^^ gcrn^) of the star as a function 
of the mass of the star. The circled point denotes the value 
at maximum mass configuration. 



that our results are in good agreement with these obser- 
vations, precisely set II and III agree with the estimates 
from PSR J0751 + 1807. In order to evaluate the radius 
of the star, we have included the BPS EoS [131 at subnu- 
clear densities which contributes to the crustal part of the 
star structure. The radius obtained for the three EoS lies 
in the range (24—26) km and the compactness parameter 
(M/R) thus lies in the range (0.05 - 0.09)Mo/fcm. From 
the observational point of view, there are large uncertain- 
ties in determination of the radius of the star [H, H^l 
which is primarily because of our lack of knowledge of the 
composition of the star atmosphere, large distance and 
also due to the presence of high magnetic fields. On the 
other hand, the general relativistic limit for the compact- 
ness parameter assuming a uniform density star with the 
causal equation of state i.e., P = e gives M/R < 4/9 [3. 
In comparison, our prediction lies far below this limit 
which reflects the softness of our EoS ably supported 
from the heavy-ion flow data (Fig. 2). Consequently, the 
flattening parameter, which is defined as the ratio of the 
polar to the equatorial radius {Rp/Re), for all the three 
case is ~ 0.59, at Qk- Rotation induces deformation in 
the shape of the star which leads to a dependence of the 
star's metric on the polar coordinate 9. In general, rota- 
tion stabilizes the star against gravitational collapse and 
therefore rotating neutron stars are more massive than 
the static ones. Further, the additional centrifugal forces 
in a rotating star help to counteract the pull of gravity, 
resulting in larger radii for a given mass. 

Moment of inertia of neutron stars plays crucial role 
in the models of radio pulsar. Independent of the rota- 
tion, i.e., slow or fast, the relation of moment of inertia 
to the matter distribution within the star is complicated. 
Among all the global properties of the star, the moment 



of inertia is the most sensitive one to the dense matter 
equation of state. In Fig. 8, we show the variation of the 
moment of inertia of the rotating neutron star as a func- 
tion of star mass for the three parameter sets. Here it 
can be seen that set I {m* /m — 0.80) results in the max- 
imum value of the moment of inertia at (145 = 7.22; 
/45 = IQ"^^ gem?) and as we move towards the higher 
effective mass values, this decrease becomes increasingly 
prominent. But for all the cases we observe that mo- 
ment of inertia of the compact stars increases rapidly as 
it approaches the Kepler velocity or till the mass shed- 
ding limit. The moment of inertia obtained for set II and 
III are 145 — 5.14 and J45 = 2.54 respectively. From 
the observational point of view, recently discovered rela- 
tivistic pulsar PSRJ0737 - 3039 feT*! could be the first 
one in which the moment of inertia is measured. Also 
the estimates from crab pulsar [b^I puts 145 > 1.61 for 
Mneb = 2.0Mq, which is the conservative estimate, and 
Ia5 > 3.04 for Mneb = 4.6M0 (latest estimate). 
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FIG. 9: Keplerian velocity (10* S*"^) as a function of the 
Gravitational mass of the star. The dot-dashed horizontal 
lines denotes the observed frequency from the pulsars. 

In Fig. 9, we plot the Kepler angular velocity (^k) 
of the star as a function of the mass of the star. The 
Kepler angular velocity obtained for the three cases are 
4141 s-\ 3999 and 3714 s"! for set I, II and III 
respectively, which corresponds to a frequency of 659IIz, 
636Hz and 591i?2:. The Kepler period, which is defined as 
Pk = 27r/J7if, obtained are 1.5277is, 1.57tos and 1.69tos 
for set I, II and III respectively. It is to be noted that 
the values that we obtained presently with the three EoS 
is in agreement with the maximum observed value of 
the Kepler period for the observed fast rotating pulsars, 
namely the PSRB1937 + 21 (i/ = 633 Hz) [63| and 
PS'i?B1957-|- 20 (i/ = 621 Hz) [6^. On the other side. 
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TABLE II: Comparison of the static (upper row) and the cor- 
responding rotational (middle row) attributes of the neutron 
star predicted with the three parameter sets of the present 
model. 



EoS 


m*/m 
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Mb 
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(Me) 


( 10"=' ffcm^) 


(ms) 


I 


0.80 


1.96 


17.5 


2.85 


2.18 










2.36 


26.4 


2.60 


2.56 


7.22 


1.52 


II 


0.85 


1.65 


16.7 


2.96 


1.81 










1.98 


25.5 


2.60 


2.11 


5.14 


1.57 


III 


0.90 


1.21 


15.0 


3.45 


1.31 










1.41 


23.7 


2.96 


1.47 


2.54 


1.69 



tabulated in Table 2. 



FIG. 10: Angular momentum of the star as a function of the 
star mass. 



the constraint imposed by general relativity on the period 
of a relativistic compact star emphasize P > 0.24 ms 
However there are indications of even faster rotating pul- 
sars from the observational point of view. Recently the 
observation of the fastest rotating neutron star, with a 
frequency;/ = 1122 i/z in an X-ray burst from the X-ray 
transient, XTE Jl 739— 285 was reported 65] in addition 
to the observation oi v = 716 Hz from J1748 — 2446ad 
[6^ , which may set interesting insight into the structural 
aspects of neutron stars. It has been mentioned that 
for the rapid rotation regime, one needs submillisecond 
pulsars with v > 1000 Hz [63] ■ However, the smaller 
Ufc, that we get has to do with the large radii that we 
obtained in the present EoS. Therefore the present EoS 
while consistent with the observed PiS'i?i?1937 -I- 21 and 
PSRB1957 + 20 data, could not account for the pulsar 
rotating with high frequency such as 1/ = 716, 1122 
This situation seems to be understandable for the fact 
that the EoS employed presently represents the class of 
softest EoS in comparison to other theoretical models 
that can explain the existence of fast rotors. In this con- 
text, EoS for skyrmionic models appear to be promising 
[6^. However, with a phase transition structure possibly 
having a crystalline color superconducting quark matter 
may account for the fast rotating pulsars |69| , before the 
mass-shedding limit. 

Similar results were obtained for the angular momen- 
tum of the star for the three cases, which is shown in 
Fig. 10. It can be seen that the angular momentum in- 
creases rapidly as the star reaches the mass-shed limit 
and finally attains a maximum value at maximum mass, 
which is shown with the filled circles for the three cases. 

Finally, for a better correlation with the underlying 
nucleon effective mass, the overall global properties for 
the static and the rotational case for the three EoS is 



V. SUMMARY AND OUTLOOK 

In our present work, we intended to present a unified 
approach to correlate the properties of dense matter EoS 
with respect to the variation in nucleon effective mass 
(m*/m = 0.80, 0.85, 0.90) defined at saturation den- 
sity and recently obtained heavy- ion collision data [40| . 
Subsequently we study the rotational attributes of the 
sequence of neutron stars in a relativistic framework in 
the mean-field approach using an effective chiral hadronic 
model generalized to include all the octet of baryons. 
The EoS employed in our present investigation satisfies 
the nuclear matter saturation properties at reasonable in- 
compressibility {K = 300 MeV). We then compared the 
resulting EoS for the three parameter sets with the ex- 
perimental flow data for both symmetric nuclear matter 
and neutron matter. We found that all the three parame- 
ter sets are in good agreement with the flow data, but on 
a critical look, parameter set II [m* /m — 0.85) satisfy 
these two combined constraint quite well, however all the 
parameter sets rather lie on the softer domain of neutron 
matter collision data. The resulting hyperon rich mat- 
ter was also found to be very sensitive to the difference 
in the nucleon effective mass. From our analysis, it was 
conclusive that hyperons form a sizeable population in 
the dense matter and with higher nucleon effective mass 
value the threshold for appearance of the respective hy- 
peron species is pushed further to the higher densities. 
Hyperons are found to have a substantial impact on the 
density region relevant to neutron star properties and is 
known to decrease the maximum mass of the star. The 
ongoing and future experimental hypernuclear programs 
such as those at Jefferson Lab, KEK, J-PARC and GSI, 
Darmstadt and FAIR etc. will provide the decisive in- 
sights into the role of these exotic forms of matter in 
dense matter EoS and hence neutron stars, so as to effi- 
ciently constrain the dense matter EoS. 

We calculated the global properties of the rotating neu- 
tron star sequences such as the mass, radius, central den- 
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sity, moment of inertia and Kepler angular velocity in 
the KEH method or the fast rotation approach. The 
maximum mass obtained at Keplerian velocity for the 
three sets lies in the range (1.4 — 2.4)Mq which is in 
very good agreement with recent observations of massive 
neutron stars such as PSR J0751 + 1807, although we 
do find substantial increase in the star radius. On over- 
all analysis, parameter set II seems to fit the observed 
mass of PSR J0751 + 1807 (2.1 ± 0.2) within la error 
bar, which can be considered as ideal parameterization in 
the present context. The compactness parameter [M / R) 
that we obtained for the three EoS represents the lower 
limit from the observational point of view. The corre- 
sponding central density at maximum mass lies in the 
range (2.6 — 3.0)po- The flattening parameter i.e., the ra- 
tio of the polar radius to the equatorial radius obtained 
for the three sets is w 0.59 which reflects the degree of de- 
formation that can result in case of rapidly rotating stars. 
Subsequently the circumferential radius obtained for the 
three cases lies in the range (23.7— 26A)Km. Further we 
obtained the Kepler period in the range (1.52 — 1.69)ms, 
which can account for the fast pulsars observed such as 
the case of PS'i?B1937-f 21 and PS'i?B1957-|-20, rotating 
with frequency of v = 633Hz and v = 62lHz respec- 
tively. However our EoS could not explain the pulsars 
rotating beyond v ~ 659Hz, keeping in mind that we 
investigated the neutron star sequences with a hyperon 
core. Possibly in order to account for pulsars rotating 
at still higher frequencies, we require a stiff EoS with 
a phase transition structure in the neutron star core as 
emphasized in Ref [69j . 

The observation of high mass stars has raised consid- 
erable amount of debate over the nature of EoS i.e., the 



stiffness or softness of the EoS at high densities. In the 
present context we would like to stress that the EoS ob- 
tained in the present model belongs to the softest pre- 
scription in comparison to other theoretical models, as 
also evident from the analysis with respect to the exper- 
imental flow constraint. We also want to emphasize that 
the present EoS is compatible with the DBHF at low 
density and the phenomenological RMF-NL3 parameter- 
ization at high densities [22[. Besides, the analysis in the 
present work throws interesting insights with regard to 
the implications of heavy-ion collision data on the neu- 
tron star structure. Here we see that parameter set I and 
HI represents the upper and the lower bound of the flow 
data in the symmetric nuclear matter case, although all 
the three sets satisfies the soft prescription of the neu- 
tron matter domain. Apparently this was also evident 
from the global properties of the star that we obtained 
in our present analysis. The present EoS results in small 
Kepler frequency before the mass-shedding limit which is 
for the fact that we obtain a large circumferential radii 
for the neutron stars in the present model, which greatly 
influences the rotational attributes of the star. However 
we would like to mention here that more systematic cal- 
culations are required before we arrive at some precise 
conclusions on the constraints of dense matter EoS and 
neutron star properties provided by the present model. 
We explored here the model of neutron star with a hy- 
peron core, it shall be interesting to study and analyze 
these results including the phase transition aspects such 
as the transition to superconducting quark matter in the 
neutron star core, which we shall take up in our next 
endevour [TOj . 
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